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ABSTRACT 

The properties of underluminous type la supernovae (SNe la) of the 91bg subclass have yet 
to be theoretically understood. Here, we take a closer look at the structure of the dim SN la 
2005bl. We infer the abundance and density profiles needed to reproduce the observed spec- 
tral evolution between — 6d and +12. 9d with respect to B maximum. Initially, we assume 
the density structure of the standard explosion model W7; then we test whether better fits 
to the observed spectra can be obtained using modified density profiles with different total 
masses and kinetic energies. Compared to normal SNe la, we find a lack of burning prod- 
ucts especially in the rapidly-expanding outer layers (v>15000km s _1 ). The zone between 
—8500 and 15000km s" 1 is dominated by oxygen and includes some amount of intermediate 
mass elements. At lower velocities, intermediate mass elements dominate. This holds down 
to the lowest zones investigated in this work. This fact, together with negligible-to-moderate 
abundances of Fe-group elements, indicates large-scale incomplete Si burning or explosive O 
burning, possibly in a detonation at low densities. Consistently with the reduced nucleosyn- 
thesis, we find hints of a kinetic energy lower than that of a canonical SN la: The spectra 
strongly favour reduced densities at > 13000km s _1 compared to W7, and are very well fitted 
using a rescaled W7 model with original mass (1.38M Q ), but a kinetic energy reduced by 
-30% (i.e. from 1.33-10 51 erg to 0.93-10 51 erg). 
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1 INTRODUCTION 

Type la supernovae (SNe la) play a key role in modern astro- 
phys ics. They are invaluable as distance indicators for cosmology 
(e.g. IPerlmutter et alj|l997l. 1 19991 : iRiess etaljfl998l : lAstier et all 
l2006l : IWood-Vasev et alj|2007h because of the high accuracy with 
which the absolute luminosity of most SNe la can be inferred. 
The luminosity varies among different objects, but the variations 
correlate with distance-independent light-curve parameters such 
as the decline in magnitudes in the B-band within 15 days af- 
ter B maximum JPhillipsI [l993h . Luminosity calibration tech- 
ni ques exploiting this fact are m ostly applied to "normal" SNe 
la ( iBranch. Fisher & Nugent! 19931) not showing poorly-understood 
peculiarities. These SNe supposedly emerge from a homogeneous 
sample of progenitors, which are thought to be C-0 white dwarfs 
(WDs) accreting matter from a non-degenerate companion star 
(single-degenerate scenario). 

In the single-degenerate paradigm, the smooth varia- 
tions among "norm al" SNe la l lBranch. Fisher & Nugentl 1 19931 : 



gration) undergoes a deflagration-detonation transition (DDT) and 
proceeds as a supersonic detonation afterwards. The efficiency and 
extent of burning in the initial deflagration may then vary from ob- 
ject to object, whi ch affects the nucleo synthesis and causes the ob- 
served variability dMazzali et al . 2007). 

Extremely sub- or superlum i nous SNe la (e.g. 



Filippenko etail 1 19921: IPhillips et ail 1 19921 : iLeibundgut et al.l 
19931 : iHowell et alj|2006t) . on the other hand, are more difficult 



lNugentetal]|l995l) can b e explained within a delayed-detonation 
scenario 1 Khokhlovl 1 1 99 lh : an initially subsonic explosion (defla- 



to explain. Here, progenitors deviating from the Chandrasekhar 
mass may play a role, or some explosions might result from a 
merger of two WDs (double-degenerate scenario). Progenitor 
systems producing peculiar SNe la might also produce some rather 
"normal" explosions, contaminating the sample of homogeneous 
explosions used for distance determination. Clarifying which 
explosion scenarios lead to SNe la at which rates is therefore 
important for supernova cosmology, but it will also be of value 
for other fields. Studies conc erned with the binary progenitor s 
and population synthesis (e.g. iRui te^ Belczynski & Fryei] l2009h . 
observed supernova rates l IGreggio. Renzini & D addi 2008) or the 
impact of supernovae on their surroundings (e.g. ISato et alj|2007l) 
will profit from understanding the origin of peculiar supernovae. 
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Thus motivated, w e analyse the 91bg-like SN 2005M 
dTaubenberger et aT 2008h. SNe of the 91bg subclass are dim and 



Filippenko et al 



Garnavich ct al 



1992l ; lLeibundgut et all 1991 ; 



20041) . They were used, with 



decline rapidly (e 
iTuratto et ai]ll99i 
other SNe, to infer the sl ope of the rel ation between luminosity 
and decline rate of SNe la dPhillipsI 19930 . but later it became clear 
that dim SNe decline even more rapidly than expecte d from a linear 
luminosity— decline-rate relatio n among normal SNe dPhillips et al .1 
ll999l : lTaubenberger et alj|2008h . Spectroscopically, 91bg-like SNe 
show characteristi c peculiarities, such as low line velocities around 
B maximum (e.g. lFilippenkoeta 71 ll992h and clear spectral s igna- 
tures of Ti II, indicating lower ionisation dMazzali et all 19971) . All 
these properties together are consistent with a low mass of newly- 
synthesised 56 Ni. To date, no elaborate explosion models have con- 
vincingly reproduced 91bg-like SNe la. Pure deflagration models 
show even lower expansion velocities than obs erved in these ob - 
jects, especially when little 56 Ni is produced ("cf. lSahu et all20o"8h . 
Delayed-detonation models mi ght explain 91bg-lik e objects within 
a unified scenario for SNe la dMazzali et al.ll2007h . Yet, there are 
hints of qualitative differen ces. One example are the improved fits 
to spectra of SN 1991bg of iMazzali etafl Jl997t> . enabled by a re- 
duction in ejecta mass and kinetic energy with respect to canonical 
values. Ultimately, only refined analyses of photometric and spec- 
troscopic properties can constrain explosion models. 

We use a spectral synthesis code to analyse the structure 
and abundance stratification of SN 2005bl, reproducing its ob- 
served spectral ev olution. The "abundance tomography" method 
dStehle et al]|2005l) . which we use, exploits the fact that the opti- 
cally thick region of the ejecta becomes smaller as time progresses. 
Thus, deeper and deeper layers contribute to spectrum formation. 
Modelling a time series of spectra, we infer the abundance profile 
from the outer envelope to as deep a layer as possible. We then 
test whether variations in mass or explosion energy are needed to 
explain the differences between spectra of normal and dim SNe 
la. This is done performing abundance tomography with various 
density profiles, and assessing the quality of the resulting spectral 
fits. The range in masses and energies sampled by the modified 
models starts at 0.5-Afch / ~5-10 50 erg and extends to 1.45-Mch / 
~2-10 51 erg (Mch' Chandrasekhar mass, IAMq). This choice has 
been motivate d by parameters inferred for observed e xtreme SNe 
la of all kinds dMazzali et al.lll997l : lHowell et alfcOOd) . 

The paper is structured as follows: First, we give a short in- 
troduction to the methods employed (Sec. [2}. We then present the 
models for SN 2005bl (Sec. [3}, discuss and assess them (Sec.Q, 
and finally draw conclusions (Sec.|5j. 



2 METHOD 

The radiative transfer code we use and the abundance tomog- 
raphy method have already been described dStehle et alj I200A 
IMazzali et al]|2008l) . Thus, we focus on aspects necessary for an 
understanding of the present study. 



2.1 Radiative transfer 

We use a ID Monte Carlo (MC) radiative transfer code 



i Abbott &Lucvll 1985 


.IMazzali & Lucvlll 9931. iLucvlll 9991. iMazzalil 


l200d and Stehle et al 


[2005J) to compute SN spectra from a given 



The code computes the radiative transfer through the SN ejecta 
above an assumed photosphere. The densities within the envelope 
are calculated from an initial density profile describing the state of 
the ejecta after homologous exp ansion has set in, which is a few 
seconds after the explosion (e.g. iRopke & Hillebranddl2005r) . The 
ejecta expand radially with r = v ■ t, where r is the distance from 
the centre, t the time from explosion (see beginning of Sec. [3}, and 
v the velocity. Radius and velocity can be used interchangeably as 
coordinates. 

From the photosphere, which is located at an adjustable v p h, 
thermal radiation [7+ = 7?„(T p h)] is assumed to be emitted 
into the atmosphere. This is of course quite a crude approxima- 
tion to the pseudo-continuous radia tion field deep in the ejecta 
dSauer. Hoffmann & Pauldrachll2006l) . Notable deviations mainly 
appear in the red and infrared, where a departure of the flux level 
from that of the observed spectra sometimes cannot be avoided. 

The radiation is simulated as "photon packets", which un- 
dergo Thomson scattering as well as line excitation-deexcitation 
processes, treated in the Sobolev approximation. The process of 
photon branching is included, which implies that the transitions for 
excitation and deexcitation can be different. In a branching event, 
the photon packet is not split up. Instead, it is emitted as a whole 
with a new frequency corresponding to a possible dow nward tran- 
sition. This "indivisible packet" approach dLucvll 19991) enforces ra- 
diative equilibrium. The downward transition is randomly selected, 
taking into account effective emission probabilities. Thus, if a large 
number of packets are simulated, the distribution of decays reflects 
the actual one. 

A modified nebular approximation, which mimics effects of 
non-local thermodynamic equilibrium (NLTE), is used to calculate 
the action of the radiation field onto the gas. For each of the 30 
zones into which the envelope is discretised here, a radiation tem- 
perature Tr and an equivalent dilution factor W are calculated. 
Thes e quantities mostly d etermine the excitation and ionisation 
state dAbbott&Lucvll 19851) . Only the variables describing the state 
of the gas are discretised; for the paths and redshifts of photons, and 
for the positions of interaction surfaces of lines, continuous values 
are allowed. 

The code iterates the radiation field and the gas conditions. 
Furthermore, T p h is automatically modified so as to match a given 
output luminosity L, taking backscattering into account. After con- 
vergence, the emerging spectrum is obtained from a formal integral 
solution of the transfer equation dLucvll 19991) . 



2.2 Density profiles 

As a first step, we adopt the density structure of the sta ndard ex- 
plosion model W7 dNomoto. Thielemann & Yokoill 1984ft as a ba- 
sis for our calculations. We then repeat the abundance tomography 
with modified density profiles, changing the total mass and kinetic 
energy of the explosion. To achieve this, the values for each grid 
point in the W7 velocity-density structure are scaled uniformly (i.e. 
all velocities by one scaling factor, and all densities by another one) 
according to: 



density and abundance profile. The aim is to infer the chemical 
structure adjusting the abundances within the envelope until an op- 
timal fit to the observed spectra is obtained. 
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Here, p' and v' are the density and velocity coordinates of each 
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Table 1. Density models used in this work, and their total kinetic energy 
E^ and mass M' . The models are named according to the scaling factors 
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grid point after the scaling. E' k and M' are the new total kinetic 
energy and mass. 

The scaled density models used in this work are listed in Ta- 
ble[T] which gives an overview of the respective and ratios. 
We have not implemented every possible energy-mass combination 
within the limits given in Section [T] Instead, we first constrained 
ourselves to a few test cases. Then, we sampled the E^—AI' plane 
more densely in the region where models of acceptable quality 
emerged (see Sec. 14. It . 

W7 naturally shows some differences with respect to more re- 
cent and realistic hydrodynamical simulations, and the scaled den- 
sity profiles can also be expected to do so. However, it is possible 
to obtain good fits to spectra of "normal" SNe la like SN 2002bo 
dStehle et al]|2005h using the W7 density structure. The results for 
the scaled profiles should therefore bring out possible differences 
between dim SNe la and normal ones. 



2.3 Abundance tomography 

The abundance tomography method dStehle et al.ll2005h uses a se- 
ries of photospheric spectra to establish the abundance distribution 
within a supernova. The idea is that the opaque core of the expand- 
ing ejecta shrinks with time. Thus, a time series of spectra carries 
information about the abundances in the ejecta at different depths. 
In the picture adopted in our code, involving an approximate pho- 
tosphere, the photosphere recedes to lower velocities with time. 
Deeper and deeper layers become visible, leaving their imprint on 
the spectra. 

The earliest spectrum available can be used to obtain the pho- 
tospheric velocity at that time and the abundances in the outer enve- 
lope. To this aim, we optimise the code input param eters to match 
that s pectrum, as in a one-zone spectral model (e.g. iMazzali et al.l 
1 19971) . The subsequent spectrum will carry the imprint of the mate- 
rial in the outer envelope and additionally that of the layers inside 
which the photosphere has receded. Because the abundances in the 
outer zone are already known, the abundances of the layers which 
have become visible can now be inferred, together with the new 
velocity of the photosphere. This procedure is then continued with 
later spectra. 

The optimum parameters inferred from a spectral model 
are subject to som e uncertainty (see also the discussion in 
IMazzali et al .120081) . One important reason for this can be degener- 



acy, which makes the spectra appear similar for different parameter 
sets. The composition adopted for an outer layer in an early-epoch 
model may therefore be in conflict with a later spectrum, if the later 
spectrum is still influenced by the outer layers. In such cases, we 
revised the parameters for the outer layers so as to optimise the 
earlier and later spectra at the same time. 



3 MODELS 

We analyse five spectra of SN 2005bl, taken at — 6d, — 5d, —3d, 
+4.8d and +12. 9d with respect to B maximum. Observational 
data and one-zone spectra l models have already been presented in 
iTaubenberger etalj | |2008h . As in that paper, we assume a total red- 
dening of E(B — V)=0.20 and a S-band rise time of 17d to cal- 
culate the time t from the onset of the explosion. 

Later spectra were not modelled, as the photosphere has 
alre ady receded to v B h< 3500km s _1 at +12. 9d [for compari- 
son, IMazzali et all d2008l) found u p h=4700km s _1 at +14d and 
w p h=2800km s" 1 at +21d in SN 2004eo]. As the photosphere 
reaches the 56 Ni-rich zone, some energy deposition should real- 
istically take place above the photosphere itself. This is is not taken 
into account in our code. Thus, to explore the innermost layers one 
would need to model nebular spectra (which are not available for 
SN 2005M) at least as a consistency check. 

The outermost ejecta of SNe la m ay partly consist of un- 



burned material (cf . Mazzali et aLl 2007 ). As the one-zone mod 



els for SN 2005b! (Taubenberger et al 



2008?) showed too much 



absorption by burned material at high velocities, we introduced a 
zone with strongly reduced abundances of burning products above 
v> 15000km s~ . This was done limiting the mass fractions of 
burning products in this zone to their value at the photosphere 
at — 6d. The unburned material at f>15000km s , which then 
constitutes >98% by mass at these velocities, is assumed to consist 
of carbon and oxygen in a ~1:1 ratio. In a preliminary stratified- 
abundance model, this was found clearly to improve the synthetic 
spectra, also with respect to the one-zone models (cf. Sec. 13. 1 ,6t . 
Consequently, we implemented such a zone in all our stratified- 
abundance models (except when using the w7e0.35 density profile, 
which has negligible densities in the outer layers). 

Below, we first discuss an abundance tomography experiment 
based on the original W7 density structure. We compare our syn- 
thetic spec tra with the observed ones and with the one-zone model 
spectra of ITaubenberger et al.l (120081) . After discussing the abun- 
dance profile, we then present models with different total mass 
and kinetic energy. Parameters (abundances, photospheric veloci- 
ties, etc.) of all models are compiled in Appendix lAl 

3.1 Abundance tomography based on W7 

The spectral models discussed here are shown in Fig.[T] where the 
most important spectral features are marked. 

3.1.1 2005 April 16: -6d, ^=8400km s" 1 

At this epoch, the supernova shows a spectrum dominated by 
singly-ionised species. In normal SNe, usuall y also doubly-ionise d 
species are detected at such early epochs (e.g. IMazzali et al .120081) . 
The zone between 8400 and 15000km s _1 is dominated by oxygen. 
The absence of the Si II A6355 emission peak suggests absorption 
by C II A6580. However, the mass fraction of C between 8400 and 
15000km s" 1 must be <10%; otherwise the C II feature would 
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Figure 1. Model se quence based on the W7 d ensity profile (05bl-w7 sequence, green lines). Observed spectra (black lines) and one-zone model spectra 
(magenta lines) fromlTaubenberger et al. 12008), which only extend to +4.8d, are shown for comparison. Identifications of the most prominent features are 
given at the beginning and the end of the sequence. Especially in the +12. 9d spectrum, many blends appear, so that the identifications are only approximate. 
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become too deep. Burned material (oxygen as a burning product 
excluded) makes up for no more than ~15% in mass according to 
the observed line depths. 

While numerous lines of intermediate-mass elements (IME) 
are visible, there are no absorptions that can unambiguously be 
attributed to Fe. We determined an upper limit to the Fe abun- 
dance of 0.01%, avoiding the appearance of a spurious Fe II fea- 
ture at ~4950A. Yet, some burning products heavier than Si and 
S are seen in the spectra: some per mille of Ti and Cr are neces- 
sary to model the absorption trough at ~4100A and the feature at 
~4700A, respectively. The se elements also c ontribute significantly 
to line blocking in the UV l lSauer et alj|2008h . 

3.1.2 2005 April 17: -5d, iv,=8100km s" 1 

The April 17 spectrum is very similar to the previous one. As the 
material directly above the photosphere is highly ionised, many fea- 
tures in this spectrum depend strongly on the abundances above 
8400km s _1 . 

At ~4950A, the stratified model has an absorption trough too 
deep. This is mostly due to the Si II A5049 line, whose strength 
largely depends on the Si abundance above i>=8400km s _1 . We 
chose this abundance so as to match the Si II A5972 line of this and 
the previous spectrum, and a simultaneous match of the Si II A5049 
line was not possible. Apart from this and some flux mismatch in 
the red, the observations are fitted well. 

3.1.3 2005 April 19: -3d, ^;,=7500km s" 1 

This model again matches the observed spectrum nicely in most 
regions. The Ti-dominated trough at ~4100A is now deeper than in 
the earlier spectra, and relatively hard to fit. A good model requires 
Ti abundances of the order of a few percent at the photosphere, and 
relatively large Ti abundances in the zones above. Thus, we set the 
Ti mass fraction to 1% between 8100 and 8400km s _1 . At larger 
velocities, the abundances are sharply constrained to some per mille 
by the features in the — 6d spectrum. 

There is still no evidence for significant amounts of Fe in the 
spectrum. Fe mass fractions of a few per cent in the layers between 
7500 and 8400km s _1 are compatible with the observations, but 
not strictly required. 

3.1.4 2005 April 26: +4.8d. iy,=6600km s _1 

In order to fit this spectrum with its low flux in the UV and blue, the 
model atmosphere must contain sufficient amounts of Ti, Cr and Fe. 
The layers at v> 8400km s _1 contain relatively small amounts of 
these elements, as dictated by the pre-maximum spectral features 
and UV flux. To compensate for this, large amounts are needed 
close to the photosphere. While the flux-blocking in the UV is quite 
sensitive to the abundances close to the photosphere, the depth of 
individual features (such as Si II A5972) is still more strongly in- 
fluenced by the composition at > 8400km s _1 . 

The most notable deviation the model from the observed spec- 
trum occurs in the blue wing of O I A7773, where there is too much 
absorption. In the outermost zone, O could only be replaced by C, 
but we already have a ~1:1 C-O mixture there. If we wanted to re- 
duce the O I absorption strength by a factor of 2 in these layers, we 
would have to postulate a ~3:1 C-0 mixture, which would seem 
quite ad-hoc. In the layers between 8400 and 15000km s _1 , the 
amount of oxygen cannot be reduced (cf. Sec. l3.1.ll . In Sec. 13.2.21 



we will show that a reduction of the density in the outer layers can 
cure this problem. 

There is some mismatch around 5700A, which seems to be 
caused by a low pseudo-continuum. This impression is however 
also due to Na I D absorption at the peak between the S II trough 
and the Si II A5972 feature. We introduced a small amount of Na 
above 8100km s _1 to obtain at least some Na I D absorption at 
+ 12. 9d. The spurious absorption appearing at +4.8d then indicates 
inaccuracies in the Na ionisation profile and its evolution wit h time, 
a common issue with synthetic spectra jMazzaliet all 19971) . 



3.1.5 2005 May 04: +12. 9d, v pfc =3250km s _1 

This model carries some conceptual uncertainty, as a possible en- 
ergy deposition by 56 Ni above the photosphere is not simulated in 
our code. Yet, the overall fit is satisfactory. Some incompatibilities 
with the abundances inferred for the outer layers could not be re- 
solved. It was, for example, impossible to get rid of the absorptions 
at ~6500 and ~7700A, which are due to Ti II AA6680, 6718, 6785 
and Si II A7849, respectively. These lines were not visible in the 
earlier spectra. 

The photosphere is now deep inside the Si-dominated zone. 
The extended red wing of the observed feature at ~9000A, caused 
mostly by Si II A9242, indicates a large Si mass fraction. On the 
other hand, the small flux in the blue and UV already demands a 
larger fraction of Fe-group elements. While the exact amounts of 
Fe, Co and Ni are somewhat uncertain, their sum can be estimated 
to be ~30%. The exact number depends on the abundances of other 
elements blocking UV flux (mostly Ti and Cr) between 3250 and 
6600km s~\ which are somewhat uncertain. 



3.1.6 Comparison to one-zone models 

Compared to one-zone models iTaubenberger et all d2008l) . shown 
in Fig. Q] as the magenta line], the stratified model sequence 
clearly constitutes an improvement in fitting the observations. The 
main reason for this is the C/O-dominated shell introduced at 
v>15000km s -1 , which makes spectral lines of burned material 
narrower. The changes with respect to the one-zone model are es- 
pecially apparent in the pre-maximum spectra: the Ca II H&K and 
Si II A6355 lines absorb less at high velocities, so that in the blue 
wings of the features only small mismatches are left. The Ti II- 
dominated trough around 4100A now has more structure. Some 
deviations, even a bit more apparent then in the one-zone models, 
remain in the red wing of Si II A6355 in the earliest spectra. This 
is largely due to re-emission in this wavelength range, caused by 
elements such as Ti and Cr which block and redistribute UV flux. 
These elements are, however, necessary to model the spectral fea- 
tures (see Sec. 13.1.11 . 



3.1.7 Abundance profile 

In Fig. [2] we compare the abundance profile deri ved in our to- 
mogr aphy experiment to the nucleosynthesis in W7 l llwamoto et al.l 
1 19991), which approximately represents a normally-luminous SN 
la jNomoto. Thielemann & Yokoi|[l984l) . In our models, unburned 
material (counting in all of the oxygen) constitutes a much larger 
fraction of the ejecta, almost the outer ~0.7Mq . Our analysis of the 
outer layers is still a bit coarse. A better-resolved analysis, yielding 
more exact results e.g. for the amount of IME between 8400 and 
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Figur e 2. Abundances of W7 nucleosynthesis calculations jlwamoto et alj 
1 19991 top panel) vs. abundance tomography of SN 2005bl, based on the 
original W7 density profile (05bl-w7 model, bottom panel). 

15000km s - , would be possible if spectra at earlier epochs were 
available (see Sec. l4.4t . 

Below the outer ~0.7M Q , the ejecta of SN 2005M are dom- 
inated by IME. The transition happens in the zone between 6400 
and 8400km s . The exact transition velocity is difficult to infer, 
as the post-maximum spectra show only a limited sensit ivity to the 
Si abundances below 8400km s '^ iMazzalietail l l 19971) have con- 
ducted a fine analysis of the O I A7773 line profile in SN 1991bg, 
and found a lower cut-off velocity of 8600km s _1 for O. We thus 
implemented a relatively sharp decrease of the O abundance in 
favour of Si below the — 6d photosphere. 

The layers between 6400 and 8400km s _1 already consist of 
~100% 56 Ni in W7. We, in contrast, find (besides IME) compar- 
atively large abundances of Ti and Cr as products of incomplete 
burning at these velocities (peak values in the order of some per 
cent). These elements contribute to the formation of the observed 
trough around ~4200A which is characteristic of 91bg-like ob- 
jects past maximum, but also to the line blocking in the UV. To 
some extent, their effects can also be mimicked by Fe, Co and Ni. 
With overly large amounts of Fe, however, individual lines in the 
optical may show up, and the flux distribution in the UV and blue 
may deviate from what is observed. Large abundances of 56 Ni or 
its decay product 56 Co outside the centre would be in conflict with 
the nebular spectra o f dim SNe la, which show very narrow lines 
dMazzali et all 19971) . 

The deepest zones that we reach with our analysis are still 



dominated by IME. However, there are signs of a transition to the 
NSE-burning zone: the large amount of line blocking and flux re- 
distribution needed to fit the +12.9d spectrum clearly points to- 
wards Fe-group abundances of several 10%. 

Compared to the one-zone models of iTaubenberger et aP 

d2008h epoch by epoch, the abundances of burning products at the 
respective photospheres are larger. In a model with homogeneous 
composition, the inferred abundances will always be some average 
between those at the photosphere and those further outwards, where 
less burning products are present. 

3.2 Models with modified density profiles 

We now show some representative spectral models based on mod- 
ified density profiles (Sec. 12.21 ). The reader interested in the abun- 
dances is referred to Section POl and Appendix [A] Here, we fo- 
cus on the differences in the spectra with respect to the W7-based 
models. To facilitate the understanding of these differences, we first 
discuss the properties of the scaled density models. 

Our scaled density models span a range of masses and kinetic 
energies (see Table [T}. Scaling the total mass and energy, the am- 
plitude and/or form of the W7 density structure is changed. What 
exactly happens depends on the Sfr ratio of the final profile with 
respect to f^J ■ Here, we distinguish the following three cases, 
which result in three classes of scaled density profiles: 

~m> = m^J- ' ^ n m ' s case ' me scaled velocity-density profile is 
obtained from W7 by reducing the density at each velocity by a 
uniform factor. The form of the density profile in velocity space is 
thus left unchanged. 

E' E 

TP < m 1 '■ Here, the energy per unit mass is reduced. This 
means that mass elements are "shifted" towards lower velocities. 
The density profile becomes steeper in velocity space, and the rel- 
ative amount of mass at high velocities is smaller. 

TP ^ m W? ' Increasing the energy per unit mass "shifts" material 
outwards, opposite to the case before. As the spectra of 91bg-like 
SNe la lack absorption at high velocities in all lines, this is gen- 
erally disfavoured. Thus, we calculated only one model sequence 
with such a density profile (w7m0.7). 

The models we discuss below are exemplary for these three 
scaling types. They are named after the underlying density models 
(e.g. 05bl-w7m0.7 is based on w7m0.7). 

3.2.1 Reduced mass and energy, -jgr = f.} V7 7 ■' 05bl-w7e0. 7m0. 7 

In these models, the density is decreased at all radii. This leads to 
a slight improvement of the spectra (Fig. [3), as the photospheres 
are deeper inside the ejecta, and the absorption velocities tend to be 
lower. Owing to the lower densities, larger mass fractions of burned 
material are necessary to fit the line depths. At high velocities, how- 
ever, oxygen still dominates and the high- velocity absorption in the 
O I A7773 line only becomes a bit weaker. 

3.2.2 Reduced energy, ^ < 05bl-w7e0.7 

In the w7e0.7 density profile (Fig. |4), the densities are sig- 
nificantly increased below ~6500km s _1 and decreased above 
~13000km s . Thus, the spectral features become narrower com- 
pared to the W7-based sequence. Line widths and positions now 
generally fit the structure of the observed spectra better. 
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Owing to the lower densities in the outer part, there is less line 
blocking by heavy elements. This decreases the flux redistribution, 
so that the flux level in the red wing of Si II A6355 and redwards is 
matched better, especially at early times. 

At the same time, the spurious high-velocity absorption is 
practically gone in Ca H&K and Si II A6355, but even more im- 
portantly in O I A7773. The reason for this is again the decreased 
density in the outer layers. As the abundances in the outer layers, 
especially of oxygen, are not fundamentally changed with the den- 
sity modification, the decrease in density translates into weaker ab- 
sorption at high velocities. 

3.2.3 Increased mass, & < ^ cx : 05bl-w7ml .25 

Despite the larger mass, the 05bl-w7ml.25 spectra show a some- 
what improved quality compared to W7 (Fig. [5] see also line ve- 
locity measurements in Sec. I4.lt . This illustrates that a super- 
Chandrasekhar total mass is not necessarily incompatible with the 
spectra of SN 2005bl. Remarkably, the improvement over the W7- 
based sequence is due to decreased densities in the outermost layers 
(«>15000km s~ ) of the warped density profile. 



3.2.4 Reduced mass, -§r > j^ 3 -: 05bl-w7m0.7 

Here, the densities in the outermost layers are increased with re- 
spect to W7. This can directly be seen in the spectra (Fig. [6j: all 
the problems which are reduced in 05bl-w7e0.7 (compared to the 
original W7-based model sequence) are now exacerbated. 



4 DISCUSSION 

4.1 Assessment of the models based on different density 
profiles - mass and kinetic energy of dim SNe la. 

Having discussed some representative cases in Section [3~2l we now 
systematically compare all models calculated on the basis of differ- 
ent density profiles. Our aim is to judge the quality of each model 
sequence in a simple and meaningful manner. To achieve this, we 
introduce three quality criteria: 

(i) Consistence of spectra. The main motivation to test modi- 
fications of the density were mismatches in the line velocities or 
widths remaining in the W7-based model sequence, especially in 
O I A7773 and Si II A6355. Other lines of the spectrum did not 
show deviations as apparent, apart from Ca II H&K, which behaves 
quite similar to Si II A635fQ. To assess if the lines are better fitted 
using different density profiles, we measured the velocities of O I 
A7773 and Si II A6355 in each synthetic and observed spectrum 
at — 6d, —3d and +4.8d. Then we calculated, for each model se- 
quence and line, the velocity difference between the observed and 
the synthetic spectra, averaged over the epochs. 

(ii) Consistence of kinetic energy. We calculated a hypotheti- 
cal kinetic energy (i?k,h yP ) for each of the abundance profiles in- 
ferred. This is the nuclear energy release (assuming a pre-explosion 



1 Ca II H&K, Si II A6355 and O I A7773 are usually the strongest lines 
in our spectra, which therefore have the highest probability of developing 
high-velocity absorptions. 



composition of equal amounts C and O) minus the binding en- 
ergy |i?bind[ of the WD (gravitational energjQ minus thermal and, 
in case of rotation, rotational energy). To judge the quality of a 
model sequence, we then compared the kinetic energy assumed in 
the density scaling (E[) to i?k,hyp. The calculation of iJk.hyp depends 
on some assumptions, the first of which is that the mass fraction 
of IME in the obscured core below the +12.9d photosphere is | 
of that above the +12.9d photosphere. Actually, this mass frac- 
tion may be between zero and the IME mass fraction above the 
+ 12. 9d photosphere. The possible error due to this is given be- 
low. The binding energy | -Ebmd | of the pr ogenitors (except for th e 
0.69Mq ones) was calculated following lYoon & Langerl ( l2005h . 
who assume a white dwarf rotation profile resulting from binary 
evolution. We used their "Bi5(Af; p c )" relation (eq. 33), assum- 
ing a central density p c of 2.0-10 9 g/cm 3 (which is typical for WD 
ignition) for M'^Mch- For sub-Chandrasekhar WDs, the central 
densities are lower even in the absence of rotation. We assumed 
negligible rotation for these cases, and ob tained the central densi ty 
for a given mass inverting formula (22) of lYoon & Langerl d2005lr[ 
(iii) Expected light-curve width. For models with good consis- 
tence based on the first two criteria, we additionally can check 
whether the density and abundance structure implies a width of the 
bolometric light curve (tlc) compatible with that of dim SNe la. We 
calculate d an expected l i ght cu rve width for each model sequence, 
following iMazzali et al.l d2007h . from the respective kinetic energy 
i?k, ejecta mass M\ and total masses of IME and NSE material 
Mime, Mnse as: 

r LC = M • k^E^M'i. 

Here, k = (O.IMime + Mnse)/ M' is proportional to the opacity 
estimate of IMazzali et al.l d2007l) . and J\f is a normalisation factor 
chosen so as to agree with their estimates of light-curve widths. In 
order to calculate k, we can assume different burning efficiencies 
in the core, as above; additionally, we may adopt as either the 
hypothetical value -Etjiyp or the value E[ from the density scaling. 
We thus calculated again an average tlc and an estimate of the 
error introduced by these degrees of freedom. In order to judge the 
models, the values tlc were compared to TLc,dim = 13.9d, which is 
the average expecte d light curve width f or the similarly dim SNe 
1991bg and 1999by dMazzali et al.l2007t) . 

We now discuss the quality of the models in terms of the three 
criteria. 

4.1.1 Line velocities — consistence of spectra 

The differences in Doppler velocity of the O I A7773 and Si II 
A6355 lines between observed and synthetic spectra are shown in 
Table|2] In this table, the models are ranked according to the abso- 
lute value of the "mean velocity difference", which is the average 
over both lines and all epochs. 

A decent match of line velocity is obtained especially for 
the 05bl-w7e0.7 model, but also, for example, for some super- 
Chandrasekhar mass models with E^/M' lower than W7. This 

2 By "binding energy" and "gravitational energy" we always mean the ab- 
solute values here, i.e. we treat them as positive numbers. 

3 This formula cannot be applied for our lowest-mass models. There- 
fore, we inferred the binding energy of a 0.69Mq progenitor from a WD 
model with constant tem perature, which uses the Timmes equation of state 
iTimmes & Arnett 1999). This equation of state takes into account a vari- 
able degree of electron degeneracy. 
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Table 2. Time-averaged line velocity differences in Si II A6355 and I A7773 from models to observed spectra (denoted by (Ai>(Si II A6355)) and 
(Au(0 IA7773)), respectively). Positive differences mean that the lines are too fast (blue) in the synthetic spectra. The models are sorted according to the 
mean velocity difference averaged among both lines (ascending in absolute value). 



Model El/E KW y M' /M w7 (Av(Si 11 A6355)) (Av(Q I A7773)) { Av) 



05bl-w7e0.7 


0.70 


1.00 


64.1 


-54.7 


4.7 


05bl-w7ml.45 


1.00 


1.45 


206.0 


86.6 


146.3 


05bl-w7eO.5m0.7 


0.50 


0.70 


-26.5 


-283.8 


-155.1 


05bl-w7e0.7ml.25 


0.70 


1.25 


-117.2 


-386.0 


-251.6 


05bl-w7e0.5m0.5 


0.50 


0.50 


292.2 


222.2 


257.2 


05bl-w7ml.25 


1.00 


125 


368.1 


359.5 


363.8 


05bl-w7eO.7m0.7 


0.70 


(1.711 


380.4 


454.9 


417.7 


05bl-w7e0.5 


0.50 


1.00 


-325.5 


-618.0 


-471.7 


05bl-w7e0.7ml.45 


0.70 


1.45 


-223.4 


-798.2 


-510.8 


05bl-w7 


1.00 


1.00 


494.7 


1012.4 


753.5 


05bl-w7el.45ml.45 


1.45 


1.45 


594.9 


942.8 


768.8 


05bl-w7e0.5ml.25 


0.50 


125 


-592.6 


-976.9 


-784.7 


05bl-w7e0.35 


0.35 


1.00 


-846.7 


-1131.7 


-989.2 


05bl-w7m0.7 


I. (10 


0.7(1 


457.9 


1815.7 


1136.8 



shows that a reduced density in the outer layers is the key to a better 
fit in the lines. To fit the observed lines well, models near the Chan- 
drasekhar mass need a E^/M' smaller by ~30— 40% with respect 
to W7. With too large a reduction in energy, line velocities become 
too low (see e.g. negative velocity differences for the 05bl-w7e0.5 
model). At low masses, generally a smaller reduction in E[/M' 
suffices: 05bl-w7e0.5m0.5 (Af'=0.69M Q ) as an extreme model 
still gives a satisfactory fit with E[/M'=(E k /M) W7 . 

Remarkably, for all mass values probed in this work, a reason- 
ably good model can be obtained (judged by the line velocities). 
The kinetic energy E' k of all well-fitting models, however, is lower 
than E^wi- 



4.1.2 Energetic consistence 

In Table [3] we show our hypothetical kinetic energy values, as well 
as the quantities from which they were calculated. We then judge 
the models by the ratio of Sk,h yp to the kinetic energy assumed in 
the density scaling (E' k ). Ideally, this ratio should be equal to one; 
the larger the deviation, the lower the rank of a model. 

For density profiles with the same mass, but different kinetic 
energy E' k , the hypothetical kinetic energy iJk.hyp usually varies sys- 
tematically. In density models with smaller E^, densities are re- 
duced in the high-velocity layers (see Sec. I3.2.2l l. which contain 
mostly unburned material. At the same time, densities are increased 
in lower layers, where the material is mostly burned. The velocity at 
which the transition (between unburned and burned material) hap- 
pens does not vary much from model to model as it is constrained 
by spectral features. Therefore, the change in the density profile re- 
sults in a larger ratio of burned to unburned material and a larger 
£k,hyp- Similarly, when E' k is increased, £k,h yP decreases. Equality, 
i.e. consistence, between -Ek.hyp and E' k is usually reached at a re- 
duced value of E' k /M' with respect to W7. The required reduction 
varies with the mass of the models (see Sec. l4.2t . 

In Table [3] two supermassive models (05bl-w7e0.7ml.25, 
05bl-w7e0.7ml.45) rank top. However, it should be noted that the 
energetic quality criterion again does not single out a certain mass, 
but sets a point of energetic consistence for each mass. All models 
with larger E^ will then feature too little nucleosynthesis to explain 
the assumed kinetic energy. The opposite holds for models with 
lower E[. 



4.1.3 Expected light-cun'e width 

We calculated estimates of the width of the bolometric light curve 
for the models ranking best in spectroscopic and energetic consis- 
tence at each mass M' . The resulting values, and those of the quan- 
tities needed for the calculation, are given in Table|4] 

The deviation of tlc from TLc,dim=13.9d strongly depends 
on the mass M'. Models with larger mass clearly tend to have a 
larger light-curve width, although they often have lower values of 
k, as relatively small abundances of burning products are needed to 
match the observed line strengths with the synthetic spectra. 

Although our expected light-curve widths are quite rough es- 
timates, one can clearly state that the criterion disfavours masses 
largely deviating from the Chandrasekhar mass. The least massive 
model, with a mass of A'/'=0.5Afch, presumably will not produce 
a broad enough light curve. Likewise, the models at A/'=1.45Afch 
will probably exhibit too broad a light curve, although these mod- 
els are not strictly incompatible with TLc.dim, as a large inaccuracy 
in 7Lc,dim results from the large mass in the obscured core. 

4.2 Location of consistent models in the E' k —M' plane 

Fig.|7]gives an overview of all models in an E^—M' plane. Accord- 
ing to their quality in spectroscopic terms, the models are marked 
with different colours; the energetic consistence is indicated by 
hatches. 

In the figure, we also indicate where spectroscopically and 
energetically consistent models can generally be expected in the 
plane: a black line is drawn approximately where the transition be- 
tween too large and too small line velocities occurs. This line is 
straight and runs from massive models with low E^/M' to sub- 
massive models with E' k /M ' k,(E^/ M) W7 . A green line approx- 
imately divides the regions of too large and too small nucleosyn- 
thetic energy yields. It lies in the same region as the line of spec- 
troscopic consistence, but is curved because the WD binding en- 
ergy shows a disproportionately strong increase with WD mass. 
For models up to M'«0.7Afch, the binding energy is negligible 
compared to the nuclear energy release, whereas at higher masses 
it is considerable, forcing the line of consistence towards smaller 
E' k and larger nucleosynthesis yields. 

The two lines of consistence are especially close to one an- 
other for masses M'<M C h- Models at M'=1.45A/ch are either 
spectroscopically or energetically inconsistent, at least under the 
assumptions we made in this work. Additionally, the light-curve 
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Table 3. Energetic balance of the models (see text for a description of the quantities). The models are ordered according to the deviation of their '/ p ratio 

fc k 

from 1 (cf. last column). 



Model 


E k/ E k,W7 


M'/M W7 


E L 


Bnucl 




^k,hyp a 


g k,hyp 

K 








[10 51 erg] 


[10 B1 erg 


[10 51 erg] 


[10 51 erg] 




05bl-w7e0.7ml.45 


0.70 


1 .45 


03)3 


2. 1 


1.15 


0.96 ± 0.10 


0.03 


05bl-w7e0.7ml.25 


0.70 


1.25 


03)3 


1.73 


0.85 


0.88 ± 0.06 


-0.05 


05bl-w7e0.5m0.7 


0.50 


0.70 


0.66 


0.87 


0.13 


0.74 ± 0.02 


0.1 1 


05bl-w7e0.7 


0.70 


1.00 


03)3 


1.26 


0.49 


0.77 ± 0.03 


-0.17 


05bl-w7e0.5m0.5 


0.50 


11.5(1 


0.66 


0.60 


0.06 


0.53 ± 0.00 


-0.20 


05bl-w7e0.5ml.25 


0.50 


1.25 


0.66 


1.89 


0.85 


1.04 ±0.07 


0.23 


05bl-w7e0.7m0.7 


0.70 


0.70 


03)3 


0.80 


0.13 


0.68 ±0.01 


-0.27 


05bl-w7ml.45 


1.00 


1.45 


1.33 


13111 


1.15 


0.75 ± 0.06 


-0.43 


05bl-w7e0.5 


0.50 


1.00 


0.66 


1.44 


0.49 


0.95 ± 0.05 


0.43 


05bl-w7ml.25 


1.00 


1.25 


1.33 


1 .56 


0.85 


0.71 ±0.04 


-0.47 


05bl-w7 


1.00 


1.00 


1.33 


1.14 


0.4') 


0.65 ± 0.02 


-0.51 


05bl-w7m0.7 


1.00 


0.70 


1.33 


0.65 


0.13 


0.53 ± 0.01 


-0.60 


05bl-w7el.45ml.45 


1.45 


1.45 


13)3 


1.69 


1.15 


0.54 ± 0.04 


-0.72 


05bl-w7e0.35 


0.35 


1.00 


0.46 


1.55 


0.4') 


1.07 ± 0.06 


1.29 



a The error estimate only reflects the error due to the unknown composition below the photosphere at +12. 9d. 



Table 4. Light-curve width estimates for the spectroscopically and energetically most consistent models at each mass M' . The models are ordered according 
to the deviation of the ratio t tlc from 1; we assume 7ix,dim = 13.9d (see text). 



Model 


Bj[/Bk,W7 


M'/M W7 


K 
[10 51 erg] 


37' 


-E-k.hyp 
[10 51 erg] 




TLc a 
Ml 


_ILC_ i 


05bl-w7e0.7 


(1.70 


1.0(1 


(13)3 


1.31 


(1.77 


0.31 ±0.10 


13.8 ± 2.5 


0.00 


05bl-w7e0.5m0.7 


0.50 


0.7(1 


0.66 


03)7 


(1.74 


0.34 ±0.10 


11.7 ± 1.9 


-0.16 


05bl-w7e0.7ml.25 


0.70 


1.25 


0.93 


1.73 


Q.88 


0.32 ±0.12 


16.1 ± 3.4 


0.16 


05bl-w7e0.5m0.5 


0.50 


(1.50 


0.66 


0.69 


(1.53 


0.47 ± 0.07 


11.2 ± 1.2 


-0.19 


05bl-w7ml.45 


1.00 


1.45 


1.33 


2.00 


(1.75 


0.28 ±0.11 


16.9 ±4.6 


0.21 


05bl-w7e0.7ml.45 


0.70 


1.45 


03)3 


2.00 


0.96 


0.35 ±0.16 


18.5 ±4.5 


0.34 



3 The error estimate reflects the errors due to the unknown composition below the photosphere at +12. 9d, and due to the uncertainties in E^. 




0.5 1.0 1.5 £'J10 51 erg] 

Figure 7. Overview of all models evaluated in this work. Colour and hatches mark the energetic and spectroscopic consistence of the models, as indicated by 

the quantities " E p h / P — 1 and (Ad), respectively (darker colours / denser hatches meaning better consistence; numerical values see Tables[2]and[3j. The green 

a k 

line divides the regions where the models have too large and too small a hypothetical kinetic energy yield, compared to the kinetic energy assumed in the 
density scaling. To the left of the black line models show too low line velocities; to the right, the opposite holds. 
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Figure 8. w7e0.7 compared to the standard W7 and the w7ml.25e0.7 den- 
sity profiles. 



criterion indicates that the width of the light curve is too large 
for these most massive models. Our least massive models with 
M'=0.5Mch are also disfavoured in this respect, as they would 
probably show too rapid a light-curve evolution. 

Criteria like those used here could give stricter limits still, if 
the chemical composition in the inner layers was known. This re- 
quires studies of nebular spectra of dim SNe la. 



4.3 05bl-w7e0.7 as a reference model 

As discussed above, our models give no clear indication for a de- 
viation from the Chandrasekhar mass. The simplest modification 
leading to better spectral fits and roughly consistent energetics is 
simply a moderate downscaling of the energy, as in the 05bl-w7e0.7 
model. Thus, we consider the 05bl-w7e0.7 model a "reference". In 
Figures [8] and [9] we show the density and abundance profiles of 
the 05bl-w7e0.7 model. Other spectroscopically consistent mod- 
els show similar densities in the outer layers, and thus also similar 
abundances in that zone. This can be verified in Figures [8] and [9] 
where the 05bl-w7e0.7ml.25 model is also plotted for comparison. 

The 05bl-w7e0.7 model features 0.46A/q of unbumed mate- 
rial (including all oxygen; C constitutes 0.04Mq). IME are dom- 
inant, with a total abundance of O.55M0 above 3350km s _1 , the 
velocity of the photosphere at +12.9d. Stable Fe is present in sig- 
nificant amounts (O.O5M0). The mass of 56 Ni (including decay 
products) above 3350km s _1 is O.O 6M0, which i s a bit higher 
than the 0.016-0.026M o found in iMazzali et ail l l 19971) above 
3500km s _1 for SN 1991bg. However, some of the 56 Ni could 
be replaced by other UV-blocking elements without changing the 
quality of the fit. Some O.23M0 of material are still hidden below 
the +12. 9d photosphere, where the IME abundances may still be 
significant (Si of the order of several 10%). 

Alternative spectroscopically consistent models show similar 
patterns in the abundance profile in velocity space, but the exact 
densities and abundances below ~10000km s^ 1 are somewhat dif- 
ferent. In 05bl-w7e0.7ml.25, as an example, the densities in the in- 
ner zones are larger. Thus, the abundances of Fe, Ti and Cr must be 
lower too keep UV opacities reasonable. For Si, moderate changes 
in the number density do not cause big changes in the spectra. 
Therefore, the smaller Fe, Ti and Cr abundances can be balanced 
by slightly larger Si abundances. 



4.4 The need for early time and nebular spectra of dim 
SNe la 

The analysis presented here could still be refined for the outer- 
most and innermost layers. The exact abundance stratification in 
the outer envelope cannot be inferred from the spectrum at — 6d, 
whose photospheric velocity is already quite low. For the inner lay- 
ers, especially the density structure and thus the abundance of Si is 
somewhat uncertain (see Sec. 14. 3b . In order to make a more precise 
study of dim SNe la possible, additional spectra in the very early 
and in the nebular phase are needed. 

The pote ntial of an analysis of the nebular spectra has already 
been shown in IMazzali et alTjl997h . Here, we would like to illus- 
trate the benefit of early time spectra, showing their sensitivity to 
the abundances in the outer envelope. We checked the influence of 
the abundances between u~11000km s _1 and 15000km s _1 on the 
— 6d spectrum, and found that these abundances have some effects 
difficult to distinguish from those of the chemical composition at 
lower velocities. Moreover, the (small) abundances of burned mate- 
rial at >15000km s _1 cannot be exactly determined, as these only 
affect the extreme blue wings of the spectral features. 

To explore the effect of the abundances in the outer envelope 
on early-time spectra, we calculated model spectra at — lOd and 

— 15d (Fig. UOt. The luminosities at these epochs were crudely 
estimated from the luminosi ty at — 6d under the assumption of a 
quadratic light curve rise (cf. lRiesseta?]|l999l) . We first calculated 
spectra assuming the 05bl-w7e0.7 density and abundance structure. 
For each of the two epochs, the photospheric position was shifted 
from its value at — 6d to higher velocities, until the backscattering 
was reasonably reduced. This resulted in photospheric velocities of 
11750 and 15200km s , respectively. 

After calculating these initial models, we explored the effect 
of changes in the chemical composition, performing three addi- 
tional code runs for each epoch. In the first two runs, we reduced 
IME and heavier elements to 20% of their original abundances, 
respectively. For the — 15d model, these changes were applied to 
the whole atmosphere. In the — lOd model, we kept the original 
composition at velocities >15200km s _1 constant, in order to show 
the sensitivity to the abundances in the zone not probed by the 

— 15d spectrum. In the third code run, finally, we removed oxy- 
gen in favour of carbon (so that the mass fraction X(C)=80%). 
This change was applied to the whole atmosphere, also at -lOd, 
as otherwise an inverted composition (larger C abundances further 
inwards) would have resulted. 

In Fig. [TO] we show the resulting spectra and give line iden- 
tifications to clarify the effect of the modified abundances. More- 
over, it is indicated which lines do and which do not change sig- 
nificantly with the modifications. At — 15d, the synthetic spectrum 
looks vastly different from the spectrum of a normal SN la. We il- 
lustrate this in the upper panel of Fig.QiJJby additionally plotting 
the earliest SN la spec trum ever observed (SN 1990N at — 14d, 
lLeibundgutetal.lll99lh . Compared to spectra of normal SNe la, 
but also to the — 6d spectrum, lines of less ionised species appear 
owing to the low temperatures, which result from the low luminos- 
ity. Si II and S II lines, which normally characterise SNe la, are 
absent. The spectrum is especially sensitive to the abundances of 
Na, Ca and Fe-group elements (with Na, uncertainties in the ioni- 
sation remain a caveat, see Sec. 13. 1 .4b . Furthermore, C I features are 
present around 6800A and 8700A. As the C abundance is already 
quite large in the outermost layers of 05bl-w7e0.7, these features 
do not react strongly to a further increase of -X'(C). However, if 
much less carbon was present, they should gradually disappear. 
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Figure 9. Abundance tomography of SN 2005bl based on w7e0.7. The abundances are plotted versus enclosed mass (panel a) and velocity (panel b). For 
comparison, we also show the 05bl-w7ml.25e0.7 abundances (panels c and d). In velocity space, the patterns of the abundance profiles are very similar. 



At — lOd, the structure of the spectrum resembles somewhat 
more that at — 6d. Yet, the spectrum has little in common with 
that of the moderately s ubluminous, spectrosc opically rather nor- 
mal SN 2004eo at - 1 Id l lPastorello et alj2007l : plotted in the lower 
panel of Fig.fTOt. Compared to — 6d, the — lOd spectrum still shows 
hints of lower temperatures: because of the scarce population of 
excited levels, the S II "W-trough" does not show up. For the same 
reason, the C II A6580 feature is weak. In the model with a larger C 
mass fraction, however, some of the strongest lines of C I begin to 
absorb at ~ 8700A. Furthermore, there are absorptions due to O I, 
Na I, Si II, Ti II, Cr II and Fe II, which should allow for an analysis 
of the abundances in the outer layers as soon as observations are 
available. 

The amount of extra information which can be inferred from 
early-time spectra will, of course, also depend on the actual lumi- 
nosity of the SN at these epochs. Larger luminosities mean higher 
temperatures, making lines of different ions appear. However, our 
results already suggest that there are interesting possibilities to in- 
fer the chemical composition of the outermost ejecta. 



5 CONCLUSIONS 

We conducted an abundance tomography of SN 2005bl and con- 
firmed that nuclear burning in dim, 91bg-like SNe la stops at less 



advanced stages compared to normal SNe la. The spectra indicate 
that the abundance of burned material above ^8500k m s -1 is much 
lower than even in moderately-luminous objects iMazzali et al.l 
l2008h . From ~8500km s" 1 down to ~3300km s"\ IME dom- 
inate the ejecta. This points towards large-scale inc omplete Si- 
burning or explosive O burning (e.g. IWoosley|[T973h . A detona- 
tion at low densitie s, as it proceeds i n the outer layers of delayed- 
detonation models dKhokhlov| [ 199 11) . may be responsible for the 
abundance pattern we find. Assuming this, we need to under- 
stand how low densities could prevail in such a large fraction of 
the envelope at the onset of the detonation. Up to now, all ex- 
plosion models wh ich pre-expand the star by a deflagration and 
then detonate (e.g. Hjllebrandt&Jjigme^ei 2000L Badenes et all 



20031, iGamezo. Khokhlov & Oranj|2004lRopke & Niemevedl 20071 
Bravo et alJl2009h jroduce larger amounts of 56 Ni. This indicates 
that either the pre-expansion is too weak or the amount of 56 Ni pro- 
duced in the deflagration stage is already too large. As it is uncer- 
tain whether a suitable single-degenerate model can be found, the 
possibility of a double-degenerate origin of dim SNe la deserves 
attention. 

Besides the abundances, we have obtained information about 
the density profile of SN 2005bl. We showed that the spectra are 
incompatible with the presence of significant amounts of oxygen 
at u>13000km s _1 . Together with the low abundances of burning 
products, this indicates a general lack of material at high veloci- 
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Figure 10. Early time synthetic spectra for 15d and lOd before B maximum, based on ()5bl-w7e0.7. The black lines are spectra calculated with the original 
composition. The other lines illustrate the most notable changes which occur when setting X(C) to 80% (red, dashed line) and when reducing IME (Mg up 
to Ca, blue, solid line) or heavier elements (Sc to Ni, magenta, dotted line) to 1/5 of their original abundances, respectively (in the velocity ranges indicated). 
For comparison, we plotted spectra of SN 1990N at -14d and of SN 2004eo at -1 Id as grey thick lines below the synthetic spectra. 

Approximate identifications are given for prominent features. Weak lines are given in parentheses; lines not reacting to changes in the abundances are shown 
with dotted marks, with the reasons for the insensitivity indicated as follows (based on a rough analysis): (b) - line is heavily blended; (o) - line is formed 
mainly at >15000km s" 1 ; (r) - relative change of abundance is too small; (s) - line is partially saturated. 



ties, a lbeit less extreme than in objects like SN 2005hk dSahu et all 
120081) . We tested whether a good fit to the observed spectra is 
possible using a modified W7 model, scaled to a different total 
mass and/or energy. Indeed, a reduction of ~30% in total kinetic 
energy yielded a spectroscopically, and also energetically consis- 
tent Chandrasekhar-mass model (05bl-w7e0.7). Such consistence 
can also be reached with somewhat super- or sub-Chandrasekhar 
mass density profiles, provided that they are similar to w7e0.7 at 
v>13000km s _1 . Deviations of >30% from the Chandrasekhar 
mass seem disfavoured. With our most massive models (1.45Mch), 
it proved impossible to obtain spectroscopic and energetic consis- 
tence at the same time. In addition, these models as well as the 
least massive ones (0.5Mch) most likely would yield a light curve 
not matching that of a dim SN la. 

Sharper constraints on density models, as well as on the abun- 
dance structure in the innermost and outermost layers may be ob- 



tained from very-early-epoch and nebular spectra of dim SNe la. 
More extensive observations are needed in order to complete our 
picture of these objects, and of SNe la in general. 
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APPENDIX A: PARAMETERS OF THE MODELS 

Table [ATI shows the code input parameters of all spectral models 
mentioned in the main paper. Apart from the abundances, the code 
takes as input the photospheric velocity v p h, the time from explo- 
sion t (see main text) and the bolometric luminosity Lboi- For differ- 
ent models of a given spectrum, these luminosities can differ a bit, 
depending on the model spectral energy distribution. In addition to 
the input values, Table I All also gives the calculated temperature of 
the photospheric black body emission, Tbb, for each model. 
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Table Al. Parameters of the models. Abundances are given only for more significant elements. 



epochs 
M 



[km s 



Element abundances ( mass fractions) 
X(C) X(0) JC(Na) X(Mg) X(AI) X(Si) X(S) X(Ca) X(Ti) X(Cr) X(Fe) a) X( 5f> 



Ni)o" 



05bl-w7e0.5m0.5 



05bl-w7e0.5m0.7 



05bl-w7e0.5ml.25 



05bl-w7e0.7m0.7 



05bl-w7e0.7ml.25 



05bl-w7e0.7ml.45 



05bl-w7m0.7 



05bl-w7ml.25 



05bl-w7ml.45 



05bl-w7el.45ml.45 



-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 
-6 
-5 
-3 
4.8 
12.9 



5.520 
5.622 
5.753 
5.870 
5.609 
5.510 
5.618 
5.736 
5.837 
5.570 
5.520 
5.618 
5.740 
5.845 
5.590 
5.525 
5.625 
5.748 
5.866 
5.602 
5.517 
5.615 
5.748 
5.861 
5.595 
5.514 
5.619 
5.744 
5.843 
5.587 
5.520 
5.617 
5.745 
5.861 
5.594 
5.523 
5.623 
5.750 
5.858 
5.594 
5.519 
5.620 
5.752 
5.864 
5.605 
5.520 
5.630 
5.745 
5.844 
5.566 
5.525 
5.627 
5.754 
5.859 
5.594 
5.515 
5.622 
5.756 
5.862 
5.598 
5.520 
5.628 
5.758 
5.861 
5.605 
5.524 
5.623 
5.760 



8200 
7975 
7600 
7050 
3350 
7600 
7275 
7000 
6400 
3250 
7900 
7600 
7100 
6550 
3350 
8300 
8050 
7550 
6900 
3350 
S250 
8000 
7650 
7250 
3350 
7950 
7600 
7100 
6550 
3250 
8350 
8100 
7600 
6X00 
3350 
S700 

8400 
7950 
7150 
3350 
8700 
8450 
8200 
7350 
3525 
8150 
7850 
7475 
6550 
3250 
8400 
8100 
7500 
6600 
3300 
8800 
8500 
7950 
7000 
3325 
8800 
8550 
8200 
7175 
3475 
9000 
8700 
8100 
7000 
3450 



9574.5 
10024.7 
10377.4 
8658.6 
10611.0 
9626.0 
10220.3 
10242.6 
8613.5 
9195.7 
9734.1 
10242.2 
10727.9 
8764.9 
9582.3 
9705.2 
10185.4 
10745.4 
8959.2 
10432.1 
10038.1 
10511.8 
10658.7 
8628.6 
10835.7 
9749.2 
10346.2 
10760.9 
8730.7 
9539.3 
9764.8 
10135.5 
10632.8 
8930.6 
10071.5 
9667.9 
10101.1 
10346.8 
8750.8 
10598.5 
9896.0 
10296.5 
10222.8 
8761.1 
10682.2 
9593.4 
9927.4 
9998.2 
8654.3 
9133.8 
9845.6 
10291.1 
10859.7 
9110.3 
9957.4 
9685.6 
10207.5 
10729.3 
8996.5 
10378.3 
9990.6 
10483.5 
10562.5 
9017.1 
10418.1 
10051.5 
10534.9 
11060.7 
9313.8 
10342.2 



0.05 
0.02 
0.00 
0.00 
0.00 
0.06 
0.03 
0.00 
0.00 
0.00 
0.06 
0.04 
0.00 
0.00 
0.00 
0.05 
0.02 
0.00 
0.00 
0.00 
0.05 
0.01 
0.00 
0.00 
0.00 
0.07 
0.02 
0.00 
0.00 
0.00 
0.06 
0.03 
0.00 
0.00 
0.00 
0.06 
0.03 
0.00 
0.00 
0.00 
0.09 
0.07 
0.00 
0.00 
0.00 
0.06 
0.03 
0.00 
0.00 
0.00 
0.05 
0.00 
0.00 
0.00 
0.00 
0.08 
0.05 
0.00 
0.00 
0.00 
0.06 
0.03 
0.00 
0.00 
0.00 
0.08 
0.09 
0.00 
0.00 
0.00 



0.87 0.0080 
0.20 0.0050 



0.15 
0.11 
0.08 
0.72 
0.17 
0.01 



0.12 
0.01 



0.87 
0.13 
0.02 



0.02 

0.01 



0.13 
0.01 



0.13 
0.01 
0.01 



0.17 
0.01 



0.81 
0.12 
0.01 



0.01 
0.01 



0.0000 
0.0000 
0.0000 
0.0080 
0.0040 
0.0000 



0.00 0.0000 

0.00 0.0000 

0.83 0.0060 

0.12 0.0030 

0.03 0.0000 

0.00 0.0000 

0.00 0.0000 

0.88 0.0070 



0.0035 
0.0000 



0.00 0.0000 
0.00 0.0000 



0.0045 
0.0045 
0.0020 



0.00 0.0000 

0.00 0.0000 

0.75 0.0050 

0.10 0.0025 



0.0000 
0.0000 



0.00 0.0000 
0.86 0.0060 



0.0030 
0.0000 



0.00 0.0000 

0.00 0.0000 

0.86 0.0045 

0.13 0.0020 

0.00 0.0000 

0.01 0.0000 

0.00 0.0000 

0.82 0.0055 



0.0030 
0.0000 
0.0000 



0.00 0.0000 

0.80 0.0060 

0.11 0.0020 

0.00 0.0000 

0.01 0.0000 

0.00 0.0000 

0.84 0.0025 



0.0013 
0.0000 



0.00 0.0000 
0.00 0.0000 



0.0040 
0.0040 
0.0000 



0.00 0.0000 

0.00 0.0000 

0.83 0.0040 

0.10 0.0020 



0.0000 
0.0000 



0.00 0.0000 
0.71 0.0043 



0.12 
0.01 
0.01 



0.0043 
0.0000 
0.0000 



(1.03 
0.03 
0.04 
0.00 
0.00 
0.13 
0.08 
0.00 
0.00 
0.00 
0.06 
0.04 
0.00 
0.00 
0.00 
0.03 
0.02 
0.00 
0.00 
0.00 
0.04 
0.04 
0.00 
0.00 
0.00 
0.10 
0.03 
0.00 
0.00 
0.00 
0.04 
0.03 
0.00 
0.00 
0.00 
0.03 
0.01 
0.00 
0.00 
0.00 
0.04 
0.04 
0.00 
0.00 
0.00 
0.05 
0.04 
0.00 
0.00 
0.00 
(1.03 
0.04 
0.00 
0.00 
0.00 
0.05 
0.05 
0.00 
0.00 
0.00 
0.06 
0.04 
0.00 
0.00 
0.00 
0.10 
0.07 
0.01 
0.00 
0.00 



0.0025 0.04 

0.0025 0.60 

0.0015 0.50 

0.0000 0.55 

0.0000 0.65 

0.0070 0.06 

0.0070 0.60 

0.0030 0.52 

0.0000 0.17 

0.0000 0.17 

0.0025 0.03 

0.0025 0.62 

0.0017 0.55 

0.0000 0.53 

0.0000 0.68 

0.0020 0.03 

0.0020 0.70 

0.0010 0.62 

0.0000 0.62 

0.0000 0.78 

0.0018 0.03 

0.0018 0.62 

0.0010 0.64 

0.0000 0.65 

0.0000 0.82 

0.0035 0.04 

0.0035 0.70 

0.0018 0.60 

0.0000 0.51 

0.0000 0.50 

0.0025 0.02 

0.0025 0.68 

0.0015 0.70 

0.0000 0.71 

0.0000 0.77 

0.0020 0.03 

0.0020 0.70 

0.0010 0.80 

0.0000 0.76 

0.0000 0.84 

0.0020 0.03 

0.0020 0.60 

0.0010 0.76 

0.0000 0.78 

0.0000 0.87 

0.0060 0.05 

0.0050 0.72 

0.0020 0.72 

0.0000 0.60 

0.0000 0.46 

0.0032 0.05 

0.0032 0.65 

0.0011 0.69 

0.0000 0.69 

0.0000 0.64 

0.0023 0.04 

0.0023 0.65 

0.0010 0.72 

0.0000 0.78 

0.0000 0.80 

0.0020 0.04 

0.0020 0.65 

0.0010 0.78 

0.0000 0.78 

0.0000 0.87 

0.0035 0.07 

0.0035 0.55 

0.0015 0.68 

0.0000 0.72 

0.0000 0.84 



0.00 
0.06 
0.13 
0.04 
0.00 
0.01 
0.01 
0.12 
0.04 
0.00 
0.01 
0.10 
0.15 
0.05 
0.00 
0.01 
0.07 
0.14 
0.06 
0.00 
0.00 
0.12 
0.16 
0.07 
0.00 
0.03 
0.07 
0.09 
0.07 
0.00 
0.01 
0.10 
0.1 I 
0.07 
0.00 
0.01 
0.10 
0.14 
0.08 
0.00 
0.01 
0.10 
0.12 
0.07 
0.00 
0.02 
0.07 
0.15 
0.05 
0.00 
0.02 
0.10 
0.12 
0.06 
0.00 
0.02 
0.08 
0.15 
0.09 
0.00 
0.01 
0.12 
0.10 
0.08 
0.00 
0.03 
0.10 
0.14 
0.07 
0.00 



0.0009 
0.0009 
0.0009 
0.0009 
0.0009 
0.0010 
0.0010 
0.0010 
0.0010 
0.0010 
0.0008 
0.0008 
0.0008 
0.0008 
0.0008 
0.0004 
0.0004 
0.0004 
0.0004 
0.0004 
0.0006 
0.0006 
0.0006 
0.0006 
0.0006 
0.0009 
0.0009 
0.0009 
0.0009 
0.0009 
0.0004 
0.0004 
0.0004 
0.0004 
0.0005 
0.0006 
0.0008 
0.0008 
0.0008 
0.0008 
0.0006 
0.0006 
0.0013 
0.0006 
0.0006 
0.0010 
0.0010 
0.0010 
0.0010 
0.0005 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0003 
0.0005 
0.0005 
0.0005 
0.0005 
0.0005 
0.0002 
0.0002 
0.0002 
0.0002 
0.0002 



0.0016 
0.0070 
0.0500 
0.0850 
0.0111 
0.0005 
0.0300 
0.1100 
0.2550 
0.2150 
0.0005 
0.0150 
0.0750 
0.1225 
0.0400 
0.0008 
0.0300 
0.0625 
0.0675 
0.0190 
0.0010 
0.0220 
0.0600 
0.0900 
0.0070 
0.0004 
0.0200 
0.0700 
0.1200 
0.0590 
0.0004 
0.0100 
0.0533 
0.0550 
0.0167 
0.0005 
0.0100 
0.0110 
0.0290 
0.0140 
0.0006 
0.0220 
0.0275 
0.0333 
0.0170 
0.0004 
0.0050 
0.0375 
0.0900 
0.0800 
0.0004 
0.0100 
0.0500 
0.0675 
0.0350 
0.0002 
0.0160 
0.0425 
0.0450 
0.0135 
0.0003 
0.0150 
0.0275 
0.0300 
0.0097 
0.0002 
0.0225 
0.0500 
0.0650 
0.0350 



0.0008 
0.0035 
0.0250 
0.0850 
0.0111 
0.0007 
0.0300 
0.1100 
0.2550 
0.2150 
0.0004 
0.0150 
0.0750 
0.1225 
0.0400 
0.0004 
0.0250 
0.0625 
0.0675 
0.0190 
0.0008 
0.0220 
0.0600 
0.0900 
0.0070 
0.0004 
0.0200 
0.0700 
0.1200 
0.0590 
0.0003 
0.0070 
0.0367 
0.0400 
0.0167 
0.0005 
0.0075 
0.0080 
0.0240 
0.0140 
0.0006 
0.0220 
0.0275 
0.0333 
0.0170 
0.0008 
0.0050 
0.0375 
0.0900 
0.0800 
0.0004 
0.0100 
0.0500 
0.0675 
0.0350 
0.0002 
0.0120 
0.0400 
0.0450 
0.0135 
0.0003 
0.0150 
0.0275 
0.0300 
0.0097 
0.0003 
0.0225 
0.0500 
0.0650 
0.0350 



0.0017 
0.0750 
0.1000 
0.1100 
0.0800 
0.0003 
0.0400 
0.1250 
0.2550 
0.1300 
0.0004 
0.0400 
0.1200 
0.1650 
0.0800 
0.0008 
0.0100 
0.1000 
0.1600 
0.0600 
0.0013 
0.0250 
0.0600 
0.0850 
0.0500 
0.0003 
0.0300 
0.1500 
0.1600 
0.1250 
0.0003 
0.0150 
0.0900 
0.1150 
0.0650 
0.0006 
0.0100 
0.0400 
0.0900 
0.0400 
0.0009 
0.0150 
0.0450 
0.0700 
0.0250 
0.0003 
0.0050 
0.0475 
0.1500 
0.1250 
0.0001 
0.0175 
0.0800 
0.1000 
0.0950 
0.0002 
0.0190 
0.0400 
0.0400 
0.0600 
0.0006 
0.0200 
0.0525 
0.0700 
0.0350 
0.0005 
0.0200 
0.0750 
0.0675 
0.0250 



0.0000 
0.0000 
0.0000 
0.0100 
0.1600 
0.0000 
0.0000 
0.0000 
0.0100 
0.2600 
0.0000 
0.0000 
0.0000 
0.0100 
0.1550 
0.0000 
0.0000 
0.0000 
0.0100 
0.1200 
0.0000 
0.0000 
0.0000 
0.0085 
0.1100 
0.0000 
0.0000 
0.0000 
0.0100 
0.2500 
0.0000 
0.0000 
0.0000 
0.0100 
0.1300 
0.0000 
0.0000 
0.0000 
0.0100 
0.0800 
0.0000 
0.0000 
0.0000 
0.0070 
0.0650 
0.0000 
0.0000 
0.0000 
0.0100 
0.2500 
0.0000 
0.0000 
0.0000 
0.0100 
0.1900 
0.0000 
0.0000 
0.0000 
0.0050 
0.1100 
0.0000 
0.0000 
0.0000 
0.0070 
0.0700 
0.0000 
0.0000 
0.0000 
0.0067 
0.0600 



a ' The abundances of Fe, Co and Ni in our models are assumed to be the sum of o6 Ni and its decay chain products (^ 6 Co and ^ 6 Fe) on the one hand, and directly synthesised / progenitor Fe on the other hand. 
Thus, they are conveniently given in terms of the l ' Ni mass fraction at t — [X ( J Ni)o], the Fe abundance at t — [X (Fe)nJ. and the time from explosion onset t. 
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